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Abstract. The neutrino energy emission rate due to for- 
mation of Cooper pairs of neutrons and protons in the 
superfluid cores of neutron stars is studied. The cases 
of singlet-state pairing with isotropic superfluid gap and 
triplet-state pairing with anisotropic gap are analysed. 
The neutrino emission due to the singlet-state pairing of 
protons is found to be greatly suppressed with respect 
to the cases of singlet- or triplet-state pairings of neu- 
trons. The neutrino emission due to pairing of neutrons is 
shown to be very important in the superfluid neutron-star 
cores with the standard neutrino luminosity and with the 
luminosity enhanced by the direct Urea process. It can 
greatly accelerate both, standard and enhanced, cooling 
of neutron stars with superfluid cores. This enables one 
to interpret the data on surface temperatures of six neu- 
tron stars, obtained by fitting the observed spectra with 
the hydrogen atmosphere models, by the standard cooling 
with moderate nucleon superfluidity. 
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1. Introduction 

Neutron star interiors contain matter of nuclear and 
supranuclear density. Various microscopic theoretical 
models predict different compositions of this matter (neu- 
trons, protons and electrons; hyperons; kaon or pion con- 
densates; quarks), different equations of state (from very 
soft to very stiff) and superfluid properties of strongly 
interacting baryonic components (nucleons, hyperons, 
quarks). 

Neutron stars are born very hot in supernova explo- 
sions but cool gradually in time. Their cooling depends 
on properties of stellar matter. Comparison of theoreti- 
cal cooling models with observational data on the surface 
temperatures of isolated neutron stars yields a potentially 
powerful method to constrain models of superdense mat- 
ter. 
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Young and middle-age (t ^ 10 4 -10 5 yr) neutron stars 
cool mainly via neutrino emission from their interiors. It 
is important, thus, to know the main neutrino produc- 
tion mechanisms. For simplicity, we restrict ourselves by 
consideration of neutron stars whose cores are composed 
of neutrons (n), with some admixture of protons (p) and 
electrons (e). The neutrino emission mechanisms in the 
stellar cores may be divided into two groups, leading to 



standard or rapid cooling (e.g., Pethick 1992). 



The standard cooling lowers the surface temperature 
to about 10 6 K in t ~ 10 4 yr. It goes mainly via the modi- 



fied Urea process (e.g., Friman & Maxwell 1979, Yakovlev 



Levenfish 1995) 
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and the nucleon-nucleon bremsstrahlung 
N + N—>N + N + v + v, 



(1) 



(2) 



where N is a nucleon (n or p) . 

Rapid cooling is strongly enhanced by the direct Urea 
reaction 



P- 



1>- 



(3) 



which can operate (Lattimer et al. 1991) only in the cen- 
tral regions of rather massive neutron stars with not too 
soft equations of state. If the direct Urea reaction is al- 
lowed, the neutrino emissivity is typically 4-5 orders of 
magnitude higher than in the standard reactions (|l|) and 
(||) , and the surface temperature decreases to several times 
of 10 5 K in t ~ 10 4 yr. 

Cooling of neutron stars can be strongly affected by 
superfluidity of neutrons and protons in the stellar cores. 
The superfluidity is generally thought to be of BCS type 
produced under nuclear attraction of nucleons. At subnu- 
clear densities p ^ po (where po = 2.8 x 10 14 g cm -3 is the 
standard nuclear-matter density) the neutron pairing oc- 
curs due to nn attraction in the 1 So state. The superfluid 
gaps depend sensitively on nn interaction model. Vari- 
ous microsc opic th eories (e.g., Tama gaki 1970| , Amundsen 
& stgaard |1985| , Baldo et al. |l992| , Takatsuka & Tama- 
gaki 1993| , 1970 ) predict these gaps to vary in the range 



2 
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from some ten keV to some MeV. However, the singlet- 
state interaction of neutrons becomes repulsive at p ~ po, 
and, therefore, the singlet-state neutron superfluidity van- 
ishes near the boundary between the neutron star core 
and the crust. Deeper in the core (p ^ po), the triplet- 
state ( 3 P2) nn interaction can be attractive to produce 
the superfiuid with an anisotropic gap. Since the number 
density of protons is much smaller than that of neutrons, 
the singlet-state pp interaction is thought to be attractive 
in the stellar core leading to proton superfluidity. The su- 
perfluidity of n and p affects the main neutrino generation 
mechanisms (0) - (||) in the neutron star cores, and hence 
cooling of neutron stars with superfiuid interiors. Super- 
fluidity always suppresses these reactions decreasing the 
neutrino luminosity of neutron stars. 

However, the appearance of neutron or proton super- 
fluid in a cooling neutron star initiates an additional spe- 
cific neutrino production mechanism associated with the 
direct interband transition of a nuclcon, 



2.1. General formalism 



N -> N + v + v. 



(4) 



The mechanism is allowed due to existence of a superfiuid 
gap in the nucleon dispersion relation. The process may be 
called the neutrino emission due to Cooper pair formation. 
It was proposed and calculated for singlet-state neutron 
superfluidity by Flowers et al. (1976). It was rediscovered 



later by Voskresensky & Senatorov (1986, 1987). Until re- 
cently, the process has been 'forgotten', and has not been 
used in the studies of the neutron star cooling. It has been 
included into recent cooling simulations by Schaab et al. 
( |1997| ) Pag e ( |1998[ ), Yakovlev et al. ( |1998| ) and Levenfish 
et al. ( |199S ) although its effect has not been analysed in 
details. 

In Sect. 2 we present the derivation of the neutrino en- 
ergy generation rate due to singlet and triplet Cooper pair- 
ing of nucleons. In the particular case of singlet-state nn 
pairing we reproduce the result by Flowers et al. ( 1976[ ). 
The cases of triplet pairing of neutrons and singlet par- 
ing of protons appear to be different. In Sect. 3 we com- 
pare the Cooper-pair neutrino emission with other neu- 
trino processes in the neutron star cores, and in Sect. 4 
we analyse the effect of the Cooper-pair neutrinos on cool- 
ing of neutron stars with superfiuid cores. 



2. Cooper-pair neutrino emissivity 

Consider neutrino emission due to Cooper pair formation 
(||) of nucleons. In the absence of superfluidity the process 
is strictly forbidden by energy-momentum conservation: 
a neutrino pair cannot be emitted by a free nuclcon. Su- 
perfluidity introduces the energy gap into the nucleon dis- 
persion relation near the Fermi surface which opens the 
reaction. 



Following Flowers et al. (1976) we will study the process 
(^) as annihilation N + N^v + voi quasinucleons N in a 
Fermi liquid with creation of a neutrino pair. The process 
goes via electroweak neutral currents; neutrinos of any 
flavors can be emitted. We will assume that nucleons are 
nonrclativistic and degenerate, and we will use the approx- 
imation of massless neutrinos. Let us derive the neutrino 
energy generation rate (emissivity) due to Cooper pairing 
of protons or neutrons in a triplet or singlet state. I n this 
way we will generalize the results by Flowers et al. ( 1976 ) 
to the case of Cooper pairing of protons, and to the most 
important case of triplet pairing. The interaction Hamil- 
tonian (h = c = 1) is given by (e.g., Friman & Maxwell 
1979b 



Gf 
2V2 



(cy Jolo — caJ-1) 



(5) 



where Gf is the Fermi weak-interaction constant, and 
the terms containing cy and ca describe contributions 
of the vector and axial vector currents, respectively. The 
factors cy and ca are determined by quark composi- 



tion of nucleons (e.g., Okun' 1990), and they are differ 



ent for n and p. For the reactions with neutrons, one 
has cy = 1 and ca — gA, while for those with protons, 
cy = 4 sin 2 0\v — 1 ~ —0.08 and ca = —9a, where 
gA ~ 1.26 is the axial-vector constant and 0\v is the 
Weinberg angle (sin 2 0w = 0.23). Notice that similar in- 
teraction Hamiltonian describes the neutrino emission due 
to pairing of hyperons in neutron star matter. The latter 
process has b een discussed in th e literature (Balberg & 
Barnea 



1998 



Schaab et al. 1 



briefly in Sect. 2.2. 

Furthermore, in Eq. (0) 

^^(1+ 7 5 )W 



and will be considered 



(6) 



is the neutrino 4-current (p— 0,1,2,3), 7^ is a Dirac matrix, 
?/>„ is a neutrino bispinor amplitude, 



(7) 



is the nucleon 4-current and <r is the Pauli vector matrix. 

The nucleon current contains 'J, the secondary- 
quantized nonrelativistic spinor wave function of quasi- 
nucleons in superfiuid matter, and <1/ + , is its Hermitian 
conjugate. ^ is determined by the Bogoliubov transfor- 
mation. The transformation for singlet-state pairing is well 



known (e.g., Lifshitz & Pitaevskii 1980). The generalized 



Bogoliubov transformation for a triplet-state 3 P2 pairing 



was studied in detail, for instance, by Tamagaki (1970) 
In the both cases \& can be written as 



* = !>* [*- lEt+ipr u av {p)< 

pa 7} 

+ e^-^V^C-p)**], 



(8) 
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where p is a quasiparticle momentum, 



E=Je 2 + A 2 p 



(9) 



is its energy with respect to the Fermi level (valid near 
the Fermi surface, at \p — pp\ <C pf), e = vp(p — pp), 
a and 77 = ±1 enumerate spin states, A p is a superfluid 
gap at the Fermi surface (A p <C p-p Vp), vp is the particle 
Fermi velocity; x<? is a basic spinor [xtXa 1 = $aa'), &pr] 
and a prj are, respectively, the annihilation and creation 
operators. U ari (p) and V ari {p) are matrix elements of the 
operators U(p) and V(p) which realize the Bogoliubov 
transformation from particle to quasiparticle states. In the 
cases of singlet and triplet pairing, one has 



U ari (p) = u p S av , IWP)| 2 = 2« P , (10) 

where 



a q 



2 V E 



(11) 



For a singlet-state pairing, the gap A p is actually indepen- 
dent of p, so that Up and v p depend only on p = \p\. For a 
triplet-state pairing, A p , u p and v p depend on orientation 
of p. Note general symmetry properties (e.g., Tamagaki 
197Cj) 



Uar,(-P) = U ar ,(p), V an (-p) = ~V va {p). 



(12) 



Let q v — (<jJv,q v ) and q' v = (ui' v ,q' v ) be 4-momenta of 
newly born neutrino and anti-neutrino, while p = {E, p) 
and p' = (E',p') be 4-momenta of annihilating quasinu- 
cleons. Using the Fermi Golden rule one can easily obtain 
the neutrino emissivity in the form 



dp dp' 



\2V2J 2" U J (2tt) 3 (2tt) 3 
dq u dq' v 



f(E)f(E') 



2u, v (2n)* 2^(2tt)3 W M + c aW,\ 
x (2n) 4 (u; v + uj' v )5^{p + p'-q u -q' u ), (13) 

where Af u =3 is the number of neutrino flavors, an overall 
factor 1/2 is introduced to avoid double counting of the 
same NN collisions, integration is meant to be carried 
out over the domain (q u + q' v ) 2 > in which the process 
is kinematically allowed, f{E) = l/[exp(E/T) + 1] is the 
Fcrmi-Dirac distribution, T is the temperature, i, k=l,2,3, 
and 



Joo = £ \(B\i> + i>\A}\ 2 , 



'in' 



(14) 



In this case \A) denotes an initial state of the quasinu- 
cleon system (the individual states {p,rj) and (p',rj') are 
occupied) and \B) is a final state of the system (the states 



(p, 77) and [p r , 77') are empty). In Eq. ( |l3| ) we neglected the 
interference terms proportional to cy ca since they vanish 
after subsequent integration over p and p'. 

Bilinear combinations of the neutrino current compo- 
nents (JgJ) are calculated in the standard manner, and in- 
tegration over dg„ and dq' v is taken with the aid of the 
Lenard integral. The result is: 



= | [ 9 V-(- 2 -q 2 )3 Q/3 ], 



(15) 



where q = (u>,q) is 4-momentum of a neutrino pair (a; = 
w v + w' v = E + E', q = q v + q' v = p + p') and g af3 is the 
metric tensor. Inserting (|l^) into ([13]) we obtain 

x {c 2 v q 2 I 00 + 4 [mklik + - q 2 ) 1] } . (16) 

where I = I xx + I yv + Izz ■ 

Now the problem reduces to 6-fold integration over 
quasinuclcon momenta p and p' within the kinematically 
allowed domain u 1 > q 2 . Since the nucleon Fermi liquid 
is assumed to be strongly degenerate, only narrow regions 
of momentum space near the nucleon Fermi-surface con- 
tribute into the reaction. Thus we can set p = pp and 
p' = pp in all smooth functions under the integral. One 
can prove that the presence of superfluidity (of energy 
gaps) makes the process kinematically allowed in a small 
region of momentum space where p is almost antiparallel 
to p' . This allows us to set p' = p in all smooth functions 
in the integrand. 

For further integration in Eq. ( [l6| ) we write dp = 
ppdkdfl and dp' = ppdk'dfl', where dfi and dfi' are 
solid angle elements, k — p — pp and k' = p' — pp. Let 
us integrate over df2' first. For this purpose, we can fix 
p and introduce a local reference frame XYZ with the 
Z-axis antiparallel to p. Let and $ be, respectively, 
azimuthal and polar angles of p' with respect to XYZ. 
Since the space allowed for p' is small, we have <C 1, 
qx — qj_COS^>, qy = q±sin^, qz = k' — k = qu, where 
q± = Pf sin sa pp 0. In this case dp' = dk' q± dq± d<I>, 



q = |p 



P'\ 2 =q 2 



q±, and 



<llo 



q± > 0, 



where qj_ = u> — qi. The quantities /oo a- n d Iik are smooth 
functions of p and p'. In these functions, we set p' = p 
which makes them independent of q. Therefore, qx and qy 
are the only variables which depend on $. The integration 
over <!> in Eq. (|l^) contains the term, 



Iik qiqk d$ = 7r 



2g|| UiUklik 



2qp zz + q± {Ixx + Iyy) 
i(5»fc - n t n k ) Iik 



(17) 



where n = p/p. Integration is performed in the local co- 
ordinate frame XYZ, but the result is transformed to the 
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dkdk' f(E)f(E') 



basic coordinate frame using tensor character of ij& . Sub- 
sequent integration over q± from to q±o is easy and 
yields 

Q {2V2J 6 (2*)* J 

X uj(uj 2 ~q\) {24(w 2 +qf l )I QQ 

+ 9v ( 5 <?f| - uj2 ) likn-iUk + 3 {uj 2 - q 2 ) / | 

Now we introduce dimcnsionlcss variables 
v-pk , vpk' E , E' 

which give 

2 7r 2 pp a/1, r 7 



Practical formulae 

Inserting @ into (|2 
physical units we have 



and returning to the standard 



157T 5 n 10 - 6 



= 1.170 x 10 21 ^ VV-TlKR 

rriN itinc 



(18) 



(19) 



erg 



(22) 



where = Pf/vf is an effective quasinucleon mass, tojv 

T/(10 9 K), k B is the Boltzmann 



is bare nucleon mass, Tg 
constant, and 



R 



1 



dtt 



z 6 dx 



(e z 



1 



(4 j 00 + 4 /) 



(23) 



2V2 

(z + z') 



6u§ 



(c z + l)(e z ' + 1) 



(2tt) 8 
(z + z') 2 



dn 



2c. 



(z + z') 2 + 



dx dx' 
(x-x') 2 

loo 



is a function to be determined. The integrand contains 
the functions ioo an( i ^ = 4i + lyy + -^zz defined by Eqs. 
@ with p' = -p (see above). From Eqs. ©-(© for 
the singlet- and triplet-state pairings we obtain the same 
expression 



<?A 



3 {z + z'f- 



, 2 ^') 2 



/n 



8w 2 u 2 



(24) 



/\2 



5 (x - x') 



likriiUk 



(20) 



V af3 ( P ) 





— « 



p 



and the integration over x and x' is restricted by the do- 
main where (z + z') > \x — x'\/vp. The outer integra- 
tion is over orientations of nucleon momentum p. Per- 
forming the inner integration over x and x' we can as- 
sume that this orientation is fixed (and the vector n is 
constant). Then the superfluid gap A p is fixed as well. 
Introducing y = A p /T we have z = \/ x 2 + y 2 and 
z' = \J x' 2 + y 2 . The integration domain can be rewrit- 
ten as (x — x 1 ) 2 — Vp(x + x') 2 < Avpy 2 /(1 — v 2 ,). In the 
nonrclativistic limit, we are interested in, vf -C 1, and the 
domain transforms to the narrow strip in the xx' plane 

near the x = x' line. It is sufficient to set x' — x and R = Cy-r s , F s = y £ 
z' = z in smooth functions and integrate over dx — x' — x 
in the narrow range \Sx\ < 2zvp- In this way we come to 
a simple equation 



In the case of the singlet-state pairing the Bogoliubov op- 
erator V(p) possesses the properties V a p(— p) = V a p(p) 
and V a p(p) = — Vf3 a (p), and has the form (e.g, Tamagaki 
[l970| ) 



(25) 



Then from Eqs. (|l4|) we have / = 0, i.e., the axial- vector 
contribution vanishes for the singlet-state pairing in ac- 
cordance with the result by Flowers et al. ( 1976 ) (to be 
exact, the main term in the nonrelativistic expansion of 
Q over (vf/c) 2 vanishes). In this case the gap is isotropic 
and integration over dfi is trivial: 



(26) 



(e* + If 



dSx 



2zvp 



(z + z') 



'\'2 



(X 



2c. 



(z + z') 



(z + z') 2 + 

(z + z'f- 
5(x- 



(x- x'f 



Fo r the triplet-state pairing, according to Tamagaki 
( |l97bT) , V af3 (-p) = -Vpafa), V a p(p) = Vp a (p), and 
V a /3 (p) = v p T a j3 (p) , where Y a p (p) is a unitary (2 x 2) ma- 
trix. Using these relationships and Eq. (10), for the triplet 
case from (llj) we obtain 



loo 



(x 



/\2 



,l\2 



hkriirik 



I = 8u 2 p ^\V a p(p)\ s 

a/3 

which yields 
R = (c 2 v + 2c 2 A )F t , 



16u 2 v 2 



(27) 



-z 5 vf (c 



v-^00 



c 2 A I), 



(21) 



Fi 



1 

4tt 



d£ly z 



z 4 dx 



(e* + l) 2 



(28) 



and the most sophisticated term containing IikniUk van- 
ishes. 



Contrary to the singlet-state paring, the axial-vector con- 
tribution docs not vanish. 
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The results @ and @ for % and 3 P 2 superfluids 
can be written in a unified manner: 



R = aF, 



(29) 



where F stands for F s or F t , while a — c v or c v + 2 c A 
is a dimensionless reaction constant that depends on the 
particle species and superfluid type. In Table 1 we list the 
values of a for singlet-state and triplet-state superfluids 
of neutrons and protons (calculated using the values of 
cv and ca cited in Sect. 2.1). (A) denotes 1 So pairing, 
while (B) and (C) are two types of 3 P2 pairing with total 
projection of the Cooper-pair momentum onto the z-axis 
equal to mj — and 2, respectively. One can hardly ex- 
pect triplet-state pairing of protons in a neutron star core 
but we present the corresponding value for completeness 
of discussion. We give also the values of a for singlet-state 
pairing of hyperons. Hyperon superflui dity h as been dis- 
cussed recently by Balberg & Barnea ( 1998 ) and incor- 
porated into c alcula tions of the neutron star cooling by 
Schaab et al. (|1998|) . The value of a for hyperons is de- 
termined by the vector constant cy of weak neutral cur- 
rents in Eq. (|^) as a sum of contributions of correspond- 
ing quarks (e.g., Okun' 1990| ). One can see that the effi- 
ciency of the neutrino emission due to singlet-state pairing 
of various particles is drastically different. The emission 
is quite open for n, S 1 * 1 , 2° but strongly (by two orders 
of magnitude) reduced for p and S~, and vanishes for A 
and E° hyperons. Notice that the values of cv, ca and a 
can be renormalized by manybody effects in dense mat- 
ter which we ignore, for simplicity. Notice also that in the 
cases of singlet-state pairing of protons and S~ the first 
non-vanishing relativistic corrections to the emissivity Q 
due to axial- vector neutral currents (~ w|) could be larger 
than the small (~ c v ) zero-order contribution of the vec- 
tor currents. Since we neglect relativistic corrections, our 
expressions for Q in these cases may be somewhat inaccu- 
rate (give reliable lower limits of Q). 



Table 1. Reaction constant a in Eq. (E9 



particles pairing 



n, E° 

n 

V, 

p 

E± 
A, E° 



'S (A) 
3 P 2 (B,C) 
'So (A) 
3 P 2 (B,C) 
'So (A) 
'So (A) 



1 + 2 5 1=4.17 
(1 - 4 sin 2 9 W ) 2 
Jw 

(2 -4 sin 2 9w) 2 




0.0064 



I -4sin 2 e w ) 2 + 2c/ 2 1 = 
= 1.17 



3.18 



Equations (|22|) and ( |26| ) for singlet-state pairing of 
neutrons with two neutrino flavors (N v — 2) were obtained 
by Flowers et al. ( 1976j ). Similar equations were derived 
by Voskresensky & Senatorov Ql986| , |1987| ). Note that the 
final expressions for Q obtained by the latter authors con- 
tain a misprint: there is it 2 instead of it 5 in the denomi- 



22 



„ 20 



CO 



18 



faC 
Sh 

^16 



14 



12 



o 1 n 14 " 3 

p = 2 xlO g cm 



T cn =l(TK 



Neutron 
pairing 




8.0 8.5 
lgT [K] 



Fig. 1. Temperature dependence of the neutrino emissiv- 
ity due to Cooper pairing of neutrons atp = 2xl0 14 g 
cm -3 and T cn = 10 9 K for superfluidity (A) (solid line), 
(B) (dashed line) and (C) (dots). 

nator although the numerical formula includes the correct 
factor 7r 5 . In addition, the expressions by Voskresensky & 
Senatorov (198£, 1987) are written for one neutrino fla- 
vor and erroneously contain the axial-vector contribution 
which is actually negligible. 

Now we obtain practical expressions for the function F 
in Eq. (|2S|). Following Levenfish & Yakovlev (1994a, b) we 
consider three types of BCS superfluid: (A), (B) and (C) as 
described above. In case (A) the superfluid gap is isotropic, 
= Aq A \t). In cases (B) and (C) the gap is anisotropic 
and depends on angle 9 between quasinucleon momen- 
tum p and the z-axis: A B = A { B) (T) y/1 + 3cos 2 6», A c = 
Aq C ^(T) sin#, respectively, where A (T) is a temperature- 
dependent amplitude. Therefore, at given T the gap 
has minimum equal to Aq (T) for quasinucleons at the 

equator of the Fermi-sphere, whereas the gap Ac has max- 

( o 

imum Aq (T) at the equator and nodes at the poles of 
the Fermi-sphere. In the cooling theories of neutron stars 
one commonly considers the nodeless pairing (B) of neu- 
trons. However, thermodynamics of nucleon superfluid is 
very model-dependent and one cannot exclude that the C- 
type superfluid appears in the neutron star cores instead 
of the B-type, at least at some temperatures and densities. 
Let us introduce the notations 



Ao(T) 
k B T ' 



T 



(30) 



where T c is the critical temperature of nucleon superfluid- 
ity. Using the standard equations of the BCS theory, Lev- 
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enfish & Yakovlev ( 1994a[ b) obtained analytic fits which 
relate Aq(T) to r at any r < 1 for all three superfluid 
types: 



va 



Vb 



= y/l-T 1.456- 



= VT^t { 0.7893 



0.157 1.764\ 

— + — J 

1.188 



vb 



(2.030- 0.4903 r 4 + 0.1727 t 8 ) 



(31) 




8.5 9.0 

IgT [K] 



Fig. 2. Temperature dependence of the neutrino emissivi- 
ties in different reactions at p = 2po for neutron superfluid 
(B) with T cn = 8 x 10 s K, and proton superfluid (A) with 



cp 



4 x 10 y K 



One can easily see that the function F in Eq. fl29| ) 
depends actually on the only parameter v and on the su- 
perfluid type. An analysis of this function from Eqs. ( p6| ) 
and (^8|) is quite similar to that carried out by Levenfish 
& Yakovlev ( |l994tj ) in their study of the effect of superflu- 
idity on the heat capacity of nucleons. Therefore we will 
omit technical details and present the final results. 

Just after the superfluidity onset when the dimension- 
less gap parameter »«1 and r = T/T c — > 1, we obtain 

F A (v) = 0.602v 2 = 5.65(1 - r), 
F B (v) = 1.204 v 2 = 4.71(1-t), 

F c (v) = 0.4013 v 2 = 4.71(1 - r). (32) 
At temperatures T much below T c one has w>l and 



Fa{v) = ^-v 



71 13/2 -2d _ 35.5 ^_ 2v 




Fig. 3. Total neutrino emissivity versus T c 
neutron star core at p = 2po and T = 10 9 K. 



i-n and T C p 



F B (v) 
F c {v) 



1.27 



4V3 
50.03 



v e 



12.1 t 2 . 



(33) 



Therefore the neutrino emission due to the nucleon pair- 
ing differs significantly from the majority of other neutrino 
reactions. The process has a threshold (becomes allowed 
at T < T c ), and the neutrino emissivity Q is a nonmono- 
tonic function of temperature. It grows rapidly with de- 
creasing T just below T c which does not happen in other 
reactions. With further decrease of T the emissivity Q 
reaches maximum and then decreases. According to Eqs. 
( |33"| ) the decrease of Q is exponential for the nodeless su- 
perfluids of (A) or (B), and it is power-law (Q cx T 9 ) for 
the superfluid (C). The power-law behaviour of Q in case 
(C) occurs due to the presence of nodes in the superfluid 
gap. At T <C T c superfluids (A), (B), and (C) suppress 
the Cooper-pair neutrino emission in the same manner in 
which they suppress the heat capacity and the direct Urea 



process (Levenfish & Yakovlev 1994a, b). 

Finally, we have calculated F(v) numerically in a wide 
range of v and fitted the results by simple expressions 
which reproduce also the asymptotes (B2) and (|33[): 



F A (v) = (0.602 v 2 + 0.5942 v 4 + 0.288 v 6 ) 



-13/2 



F B (v) 



x (^0.5547 + v/(0.4453) 2 + 0.01130 v 2 

x exp (-y/Av 2 + (2.245) 2 + 2.245 

1.204 v 2 + 3.733 w 4 + 0.3191 v 6 
1 + 0.3511 v 2 



1/2 
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(o.7591 + V(0-2409) 2 + 0.3145 v- 



x exp {-^jAv 2 + (0.4616) 2 + 0.4616 

F c (v) = (0.4013 v 2 - 0.043 v 4 + 0.002172 v e ) 

x (1- 0.2018 v 2 + 0.02601 v 4 

- 0.001477 w 6 + 0.0000434 v 8 )^ 1 . 



(34) 



The maximum fit error is about 1% at v ~ 4 for Fa(v); 
about 3.4% at v 2 for Fb(v); and about 3% at » « 1 
for i*b(i>). 

3. Efficiency of Cooper-pair neutrinos 

Equations (p2|), ( p9| ) and (|34|) enable one to calculate the 
neutrino emissivity Q due to Cooper pairing of nucleons or 
hyperons. For illustration, we use a moderately stiff equa- 
tion of state in a neutron star core proposed by Prakash 
et al. ( |1988|) (the version with the compression modulus 
Kq = 180 MeV, and with the same simplified symme- 
try f actor Sv that has been used by Page & Applegate, 
1992] ). According to this equation of state, dense matter 



consists of neutrons with admixture of protons and elec- 
trons (no hyperons). We set the effective nucleon masses 
m* N = 0.7 tun , for simplicity. Since the critical tempera- 
tures T cn and T cp of the neutron and proton superfluids 



are model dependent (e.g., Tamagaki 197C, Amundsen & 
Cstgaard |1985| , Baldo et al. |1992| , Takatsuka & Tamagaki 
1992, 1970D , we do not use any specific microscopic super- 



fluid model but treat T cn and T cp as free parameters. 

Figure |l| shows temperature dependence of the emis- 
sivity Q produced by Cooper pairing of neutrons in the 
neutron star core at p = 2 x 10 14 g cm -3 . The critical tem- 
perature is assumed to be T cn = 10 9 K. Given p is typical 
for the transition from the singlet-state to the triplet-state 
pairing (Sect. 1). Thus various superfluid types are possi- 
ble according to different microscopic theories. We present 
the curves for all three superfluid types (A), (B) and (C) 
discussed in Sect. 2.2. A growth of the emissivity with 
decreasing T below T c is very steep. The main neutrino 
emission occurs in the temperature range 0.4 T c T ^ T c . 
In this range, the emissivity depends weakly on the super- 
fluid type and is rather high, of the order of or even higher 
than in the modified Urea reaction ([!]) in non-superfluid 
matter. This indicates that Cooper-pair neutrinos are im- 
portant for neutron star cooling (Sect. 4). 

Figure || shows temperature dependence of various 
neutrino energy generation rates in a neutron star core 
at p = 2pQ. The neutron critical temperature is assumed 
to be T cn — 8 x 10 8 K, while the proton critical tem- 
perature is T cp = 4 x 10 9 K. We show the contributions 
from the modified Urea reaction ([!]) (sum of the proton 
and neutron branches), nucleon- nucleon bremsstrahlung 
(^) (sum of nn, np and pp contributions), Cooper pair- 
ing of nucleons (sum of n and p reactions), and also the 
total neutrino emissivity (solid line) . The direct Urea pro- 
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Fig. 4. Temperature dependence of the neutrino emissiv- 
ities in main reactions at p = bpo for neutron superfluid 
(B) with T cn — 8 x 10 s K, and proton superfluid (A) with 



T, 



Cp 



4 x 10 y K 



cess (||) is forbidden at given p. Here and in what follows 
the rates of the reactions (p])-(p[ ) are taken as described 



in Levenfish & Yakovlev (1996), with proper account of 



suppression of the reactions by ne utron a nd/or proton su- 
pcrflui dity ( Levenfish & Yakovlev 1994a , Yakovlev & Lev- 
enfish 1995 ). The neutron superfluid is assumed to be of 
type (B) and the proton superfluid of type (A). A large 
bump of the total emissivity at T ~ 10 8 8 K is produced 
by the Cooper pairing of neutrons. If neutrino emission 
due to this pairing were absent the total neutrino emissiv- 
ity at T <, T cn would be 2 - 4 orders of magnitude smaller 
owing to the strong suppression by the nucleon superflu- 
idity. The Cooper pairing can easily turn suppression into 
enhancement. 

Figure 3 shows a joint effect of the neutron and proton 
superfluids onto the total neutrino emissivity Q at p = 2po 
and T = 10 9 K. Presentation of Q as a function of T cn and 
T cp allows us to display all the effects of superfluidity onto 
the neutrino emission. If T > T cn and T > T cp we have 
Q in nonsuperfluid matter (a plateau at small T cn and 
T cp ). For other T, the neutron and/or proton superfluid- 
ity affects the neutrino emission. It is seen that the neu- 
tron pairing at T ^ T cn may enhance the neutrino energy 
losses. A similar effect of the proton pairing at T ^ T cp 
is much weaker as explained above. In a strongly super- 
fluid matter (highest T cn and T cp ) the neutrino emission 
is drastically suppressed by the superfluidity. 

Figure ^| shows temperature dependence of some par- 
tial and total neutrino emissivities for the same T rn and 
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Fig. 5. Domains of T cp and T cn , where different neutrino emission mechanisms dominate (different shaded regions). 
In addition to the mechanisms we include also the neutrino emission due to ee bremsstrahlung (Kaminker et 

al. 1997). Figure (a) corresponds to the standard neutrino emission at p = 2 po and T — 3 x 10 s K, while Figure (b) 
refers to the emission enhanced by the direct Urea process at p = 5 po and T = 10 s K. Dashed lines correspond to 
T cn = Tor T cp = T. 



T cp as in Fig. g but for higher density p — 5po- The equa- 
tion of state we adopt opens the powerful direct Urea pro- 
cess at p > p cr = 1.30 x 10 15 g cm~ 3 = 4.64 p . We 
take into account all major neutrino generation reactions 
(0) -(|])- For simplicity, we do not show all partial emis- 
sivities since the total emissivity is mainly determined by 
the interplay between the direct Urea and Cooper pairing 
reactions. 

The effect of Cooper-pair neutrinos is seen to be less 
important than for the standard neutrino reactions in Fig. 
H but, nevertheless quite noticeable. It is especially pro- 
nounced if T ^ T cn <§; T cp . In this case, the strong proton 
superfluid suppresses greatly the direct Urea process, and 
the emission due to Cooper pairing of neutrons can be 
significant. 

Finally notice that while calculating neutron star cool- 
ing one often assumes the existence of one dominant neu- 
trino generation mechanism, for instance, the direct Urea 
process for the enhanced cooling or the modified Urea pro- 
cess for the standard cooling. This is certainly true for 
non-superfluid neutron-star cores, but wrong in super- 
fluid matter. In the latter case, different mechanisms can 
dominate (Fig. ||) at various cooling stages depending on 
temperature, T cn , T cpi and density. 

In particular, the Cooper-pair neutrinos dominate the 
standard neutrino energy losses at T ^ 10 9 K for a mod- 
erate neutron superfluidity (0.12 < T/T cn < 0.96). This 
parameter range is important for cooling theories. At the 



early cooling stages, when T <; 10 9 K, the Cooper-pair 
neutrinos can also be important but in a narrower range 
around T/T cn fa 0.4 especially in the presence of the pro- 
ton superfluid. As mentioned above, Cooper-pair neutri- 
nos can dominate also in the regime of rapid neutrino 
emission if the nucleons of one species are strongly su- 
perfluid but the other nucleons are moderately or weakly 
superfluid (see Fig. ||d). 

Figures ||a and b display the domains of T cn and T cp , 
where different neutrino mechanisms are dominant. Fig- 
ure ||a refers to the standard cooling at p — 2 po and 
T = 3 x 10 s K. Figure corresponds to the rapid cool- 
ing at p — 5po and T = 10 8 K. Dashes show the lines of 
T = T cn and T = T cp , which separate T cn -T cp planes into 
four regions. In the left down corners enclosed by these 
lines, matter is nonsuperfiuid. The right down corners cor- 
respond to superfluidity of n alone; the left upper corners 
to superfluidity of p alone, and the right upper corners 
to joint n and p superfluidity. One can observe a variety 
of dominant cooling mechanisms regulated by the super- 
fluidity. If both n and p superfluids are very strong, they 
switch off all the neutrino emission mechanisms involving 
nucleons (and discussed here). In this regime, a slow neu- 
trino emission due to ee-bremsstrahlung (Kaminker et al. 



1997) survives and dominates. This mechanism is rather 



insensitive to the superfluid state of dense matter. 
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lg t [yr] lg t [yr] 

Fig. 6. Redshifted surface temperature T s °° versus age t for the standard (a) and enhanced (b) cooling of the 1.3 M© 
and 1.48 M© neutron stars, respectively. Dotted curves are for non-superfluid stars. Solid and dashed curves are for 
superfluid stars (marked with (lgT c „, \gT cp )) including and neglecting Cooper-pair neutrinos, respectively. Solid and 
dashed curves (7.8,7.8) in Fig. (a), and (8.3,8.6) and (10,8.4) in Fig. (b) coincide. 



1988D as has Yakovlev 



4. Models of cooling neutron stars 

To illustrate the results of Sects. 2 and 3 we have per- 
formed simulations of neutron-star cooling. We have 
used the s ame c ooling code as described by Levenfish & 
Yakovlev (|l996|) . The general relativistic effects are in- 
cluded explicitly. The stellar cores are ass umed to have 
the same equation of state (Prakash et al 
been used in Sect. 3. The maximum neutron-star mass, 
for this equation of state, is 1.73 M©. We consider the stel- 
lar models with two masses. In the first case, the mass is 
M = 1.3 Mq, the radius R = 11.87 km, and the central 
density p c — 1.07 x 10 15 g cm" 3 is below the threshold 
(per) of the direct Urea process; this is an example of 
the standard cooling. In the second case M = 1.48 M©, 
R = 11.44 km, and p c = 1.376xl0 15 g cm" 3 . The powerful 
direct Urea process is open in a small central stellar ker- 
nel of radius 2.32 km and mass 0.035 M©, producing en- 
hanced cooling. Notice, that in calculations of the equation 
of state of the stellar core in our earlier articles (Levenfish 
& Yakovlev 1996] , and references therein) the parameter 
no (standard saturation number density of baryons) was 
set equal to 0.165 fm" 3 . In the present article, we adopt 
a more natural choice rio = 0.16 fm~ 3 and use the model 
of the rapidly cooling star with somewhat higher mass 
than before (1.44 M©). The nucleon effective masses are 
set equal to 0.7 of their bare masses. 

For simplicity, the nucleons are assumed to be super- 
fluid everywhere in the stellar core. We suppose that the 



proton superfluidity is of type (A), while the neutron su- 
perfluidity is of type (B). We make the simplified assump- 
tion that T cn and T cp are constant over the stellar core 
and can be treated as free parameters (see Sect. 3). 

Our cooling code includes the main traditional neu- 
trino reactions in the neutron star core (|l|)-(||), suppressed 
properly by neutron and proton superfluids (Levenfish & 
1996), supplemented by the Cooper neutrino 



emission by neutrons and protons (Sects. 2 and 3). In ad- 
dition we include the neutrino emission due to electron- 
ion bremsstrahlung in the neutron star crust using an ap- 



proximate formula by Maxwell (1979). The neutron-star 
heat capacity is assumed to be the sum of the capaci- 
ties of n, p, and e in the stellar core affected by n and 
p superfluids (Levenfish & Yakovlev 1994b| ). We have ne- 
glected the heat capacity of the crust due to small crustal 
masses for the chosen stellar models. Our cooling code is 
based on approximation of isothermal stellar interior valid 
for a star of age t > (10-10 3 ) yr, inside which the ther- 
mal relaxation is over. We use the relationship between 
the surface and interior stellar temperature obtained by 
Potekhin et al. ( 1997 ). We assume that the stellar atmo- 
sphere may contain light elements. Then we can compare 
our results with observations of thermal radiation inter- 
preted using the hydrogen or helium atmosphere models. 
However the mass of light elements is postulated to be 
insufficient (^S 10 _14 Af©) to affect the cooling. 

Figures ||a and b show typical cooling curves (depen- 
dence of the effective surface temperature T s °° versus stel- 
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Fig. 7. Standard cooling curves (1.3 M ) compared with 
observations (Table ^). Error bars are 90%-95% estimates 
of T s °° obtained by fitting the observed radiation spectra 
with black-body spectrum (filled circles) and atmosphere 
models (open circles) . Dotted curve corresponds to a non- 
superfluid star. Solid and dashed curves labelled as in Fig. 
H are for superlfuid stars. Dashed region is formed by stan- 
dard the curves for different T cn and T cp . 

lar age t) for the neutron stars with superfluid cores. 
The effective temperature is redshifted (as detected by 
a distant observer). Figure ||a displays the standard cool- 
ing, while ^3 shows the enhanced cooling. It is seen that 
Cooper-pair neutrino emission can be either important 
or unimportant for the cooling of both types depending 
on T cn , T cp and t. Its effect is to accelerate the cooling, 
and the stronger effect takes place if the neutron or pro- 
ton superfluidity is switched on in the neutrino cooling 
era (t ^ 10 5 yr). As mentioned in Sect. 3, the neutrinos 
provided by neutron pairing can dominate the neutrino 
emission from the modified Urea process at temperatures 
T T cn < 10 9 K. Therefore, if the neutron superfluidity 
with T cn 10 9 K is switched on in the stellar core, the 
standard cooling is accelerated. The acceleration is espe- 
cially dramatic if the modified Urea was suppressed by 
the proton superfluidity before the onset of the neutron 
superfluidity (T cn <C T cp ). Then the slow cooling looks 
like the enhanced one (cf Figs, ^a and b). On the other 
hand, Cooper-pair neutrinos can be unimportant. The ex- 
ample is given by the curve (lgT c „ = 7.8, lgT cp = 7.8) in 
Fig. ^a. In this case, the star enters the photon cooling era 
(with the photon surface luminosity much larger than the 
neutrino one) with the internal temperature T ^ 10 8 K. 
The superfluidity appears later, and has naturally almost 



no effect on the cooling. If T cn <C T cp the Cooper neu- 
trino emission by neutrons can dominate even the powerful 
direct Urea process and accelerate the enhanced cooling 
(Fig. ^b). In other cases (e.g., the curves (lgT c „ = 8.3, 
lgT cp = 8.6) and (lgT cn = 10, lgT cp = 8.4)) the effect 
of Cooper-pair neutrinos on the enhanced cooling can be 
unimportant. 

Figure [7] compares the standard cooling curves 
(1.3 Mq) with the available observations of thermal radi- 
ation from isolated neutron stars. The observational data 
are summarized in Table 2. We include the data on four 
pulsars (Vela, Geminga, PSR 0656+14, PSR 1055-52) and 
three radioquiet objects (RX J0822-43, IE 1207-52, RX 
J0002+62) in supernova remnants. The pulsar ages are 
the dynamical ages except for Vela, where new timing 
results by Lyne et al. (1996) are used. Ages of radio- 



quiet objects arc associated to ages of their supernova 
remnants. The error bars give the confidence intervals of 
the redshifted effective surface temperatures obtained by 
two different methods. The first method consists in fitting 
the observed spectra by neutron-star hydrogen and/or he- 
lium atmosphere models (open circles), and the second one 
consists in fitting by the blackbody spectrum (filled cir- 
cles). Dashed region encloses all standard cooling curves 
for a 1.30 M star with different T cn and T cp (from 10 6 K 
to 10 10 K). Notice that the Cooper-pair neutrinos intro- 
duce into the cooling theory some new "degree of free- 
dom" which helps one to fit the observational data. For 
instance, we present a solid curve which hits five error 
bars for the atmosphere models at once. We would not 
be able to find a similar cooling curve if we neglected the 
Cooper-pair neutrino emission. The dashed curve in Fig. [?] 
shows that Cooper-pair neutrinos make the standard cool- 
ing very sensitive to the superfluid parameters. A minor 
change even of the one superfluid parameter yields quite 
a different cooling curve, which hits two error bars only. 
Such a sensitivity is important for constraining T cn and 
T cp from observations. 

Thus, the majority of observations of thermal radia- 
tion from isolated neutron stars can be interpreted by the 
standard cooling with quite a moderate nucleon super- 
fluidity in the stellar core. These moderate critical tem- 
peratures do not contradict to a wealth of microscopic 
calculations of T cp and T cn . Notice that it is easier for us 
to explain the "atmospheric" surface temperatures than 
the blackbody ones. This statement can be considered as 
an indirect argument in favour of the atmospheric inter- 
pretation of the thermal radiation. Although the theory of 
neutron-star atmospheres is not yet complete (e.g., Pavlov 



& Zavlin 1998) the atmospheric fits give more reasonable 
neutron-star parameters (radii, magnetic fields, distances, 
etc.) which are in better agreement with the parameters 
obtained by other independent methods. 

Even in our simplified model (one equation of state, 
two fixed neutron-star masses, constant T cn and T cp over 
the stellar core) it is possible to choose quite definite su- 
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Table 2. Observational data 



Source 


igt, Est] 


Atmosphere model b ' 


l mOO Vt s~l 

lgT s , [K] 


Confidence level ' 


Reference 


RX J0822-43 


3.57 


Hydrogen atmosphere 


c> 90+O.O2 


95.5% 


Zavlin et al. ( 


1998b 


) 


Black body 


a ci +0.05 
6.6II0.O6 


95.5% 


Zavlin et al. ( 


1998b 


) 


IE 1207-52 


3.85 


Hydrogen atmosphere 


n 1 ,-,+0.05 

6-lOIo.oe 


90% 


Zavlin et al. ( 


1998a 


) 


Black body 




90% 


Zavlin et al. ( 


1998a 


) 


RXJ0002+62 


3.95 c) 


Hydrogen atmosphere 


a no+0.03 

6-03I . 03 


95.5% 


Zavlin et al. ( 


1998b 


) 


Black body 


n 1 o 4-0 18 

6-18I0.I8 


95.5% 


Zavlin et al. ( 


1998b 


) 


PSR 0833-45 
(Vela) 


4 4 d ) 


Hydrogen atmosphere 


t- nn+0.04 


907o 


Page et al. (1996) 


Black body 


6 94+ 03 




Ogelman ( 


1995 


) 


PSR 0656+14 


5.00 


Hydrogen atmosphere 


r 79 +0.04 
J -' -0.02 




Anderson et al. (1993) 


Black body 


5 96+ 02 

J - yu -0.03 


90% 


Possenti et al. 


1996) 


PSR 0630+178 
(Geminga) 


5.53 


Hydrogen atmosphere 


c 9 r+0.08 


90% 


Meyer et al. ( 


1994) 


Black body 


r 7 r+0.05 
°- ' d -0.08 


90% 


Halpern, Wang (1997) 


PSR 1055-52 


5.73 


Black body 






Ogelman ( 


1995 


) 



Confidence level of T s °° (90% and 95.5% correspond to 1.64cr and 2a, respectively); dash means that the level is not 
indicated in cited references. 
b ' Method for interpretation of observation. 



^ The mean age taken acco rding, to Craig et al. (1997) 
d) According to Lyne et al. (1996). 



perfluid parameters to explain most of observations. How- 
ever, one needs more elaborated models of cooling neutron 
stars to obtain reliable information on superfluid state of 
the neutron star cores. We expect to develop such models 
in the future making use of the results obtained in the 
present article. 



5. Conclusions 

We have reconsidered the neutrino emission rate due to 
Cooper pairing of nucleons in the neutron star cores (Sect. 
2). We have presented the results in the form convenient 
for practical implications in three cases: (A) singlet-state 
pairing, (B) triplet-state 3 P2 pairing with zero pro- 
jection of the Cooper pair momentum (mj = 0) onto 
quantization axis, and (C) triplet-state pairing with max- 
imum (mj — 2) momentum projection. For the singlet- 
state pairing o f neu trons, our results agree with those by 
Flowers et al. ( 1976| ). For the triplet-state pairing our con- 
sideration is original. Notice an essential difference of the 
Cooper-pair neutrino emissivities for singlet and triplet- 
state super-fluids and also for neutrons and protons. In 
Sect. 3 we have analysed the efficiency of the Cooper-pair 
neutrino emission at different densities in the neutron star 
cores as compared with the traditional neutrino produc- 
tion mechanisms including a powerful direct Urea process 
allowed at high densities. Contrary to the non-superfluid 
cores where the main neutrino emission is produced either 
by the modified or by the direct Urea processes (depending 



on equation of state and density), very different neutrino 
mechanisms can dominate in the superfluid cores at cer- 
tain temperatures T and superfluid critical temperatures 
T cn and T cp . In particular, neutrinos produced by pair- 
ing of neutrons can be very important if T ^ T cn -C T cp . 
The importance of these neutrinos in the standard and 
rapid cooling of the neutron stars has been analysed in 
Sect. 4. We show that, under certain conditions, neutri- 
nos provided by pairing of neutrons can greatly accelerate 
both standard and enhanced cooling of middle-age neu- 
tron stars (t= 10 4 -10 5 yr). In particular, the accelerated 
standard cooling can mimic rapid cooling of the stars. The 
Cooper-pair neutrinos modify the cooling curves and en- 
able us to explain observations of thermal radiation of 
several neutron stars by one cooling curve at o nce. T his 
confirms the potential ability (Page & Applegate, 1992|) to 
constrain the fundamental parameters of superdense mat- 
ter in the neutron star cores, the critical temperatures of 
neutron and proton superfluids, by comparing theory and 
observation of neutron stars. 
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